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Abstract 

We present expressions in terms of a double infinite series for the Stieltjes 
constants 7Jt(a). These constants appear in the regular part of the Laurent 
expansion for the Hurwitz zeta function. We show that the case 7fc(l) = 7 
corresponds to a series representation for the Riemann zeta function given much 
earlier by Brun. As a byproduct, we obtain a parameterized double series 
representation of the Hurwitz zeta function. 
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Introduction and statement of results 

Herein we supplement our very recent work [5] showing how a method of Addison 
[2] for series representations of the Euler constant 7 may be generalized in many 
different ways. We focus on the Stieltjes constants 7fc(a) and recall the defining 
Laurent expansion [HI El El [16] 

C(s,a) = _ + £L_J_ 7n ( a )( a _i)» s ^l. (1.1) 

where £(s, a) is the Hurwitz zeta function. For Re s > 1 and Re a > we have 

00 j 

C(s,a) = V- -, (1.2) 

and by analytic continuation ((s, a) extends to a meromorphic function through out 
the whole complex plane. Notationally, we let ((s) = ((s, 1) be the Riemann zeta 
function [9J [HI [131 US] and ip = T'/T be the digamma function (e.g., pfl [31 [10]) with 
the Euler constant 7 = 7o(l) = — V'(l)- In addition, Pi(x) = Bi(x — [x]) = x—[x] — l/2 
denotes the first periodized Bernoulli polynomial, with {x} — x — [x] the fractional 
part of x. 

The sequence {7fc(a) }£L exhibits complicated changes in sign with k. For in- 
stance, for both even and odd index, there are infinitely many positive and negative 
values. Furthermore, there is sign variation with the parameter a. These features, 
as well as the exponential growth in magnitude in k, are now well captured in an 
asymptotic expression ([12], Section 2). In fact, though initially derived for large 
values of k, this expression is useful for computational approximation even for small 
values of k. 
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In this paper, we give exact expressions for the Stieltjes constants as double infinite 
series, and one of the series has an exponentially fast rate of convergence. In a 
particular case for a = 1, we are able to recover a result that corresponds to an earlier 
series representation for the Riemann zeta function given by Brun [4j. Brun performed 
'horizontal' and 'vertical' summations, and so our approach is very different. 

Proposition 1. Let Re a > 0. Then we have 

7o (a) = -^(a) = - In a+- + EE 73^, ^ 

u n =l j=l J 

and for t > 1 



- + -Wa + EE iZ ^^fa + A). (L4 ) 



As a byproduct of our proof of Proposition 1, we obtain the following double series 
representation. We have 

Corollary 1. For Re s > 1 and Re a > 0, we have 

C(s, a) = a~ s + + Y 2 n ( s - 1) y , 1 J x . (1.5) 

Brun [1] realized the following representation of the Riemann zeta function for Re 
s > 1: 

C( s ) = J_ + l_/3( s ), (1.6a) 
s — 1 

where 

(L66) 

Therefore, we recover his result as the special case at a = 1. 



Corollary 2. We have 
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-| oo 

-Y 

* n=0 



in+1 



- ^(2 n ) 



= $:^(2" +1 )-^(2")-ln2]. 

n=0 

Corollary 3. We have for Re a > 

oo 

7o (a) = - In a + £ [</>(a2 n+1 ) - ^( a 2«) 



(1.7) 



In 2]. 



(1.8) 



n=0 



This results agrees with (2.24) of [5]. 

Proposition 2. We have for I > 0, Re a > 0, and integers k > 2 



lh/a h/ +1 



1 /l 



2 a 



£ J,n El 2 U 1 



n=0 ^ j=0 



ln'(6j + a) h/(6(j + 1) + a 
(fej + a) + 



k-l 



+ E 



h/[6(j + m/fc) + a] 



fc^i [6(j +m/fc) + a] 



(b(j + 1) + a) 

(1.9a) 



ln £+1 a 



oo fe— 1 



+ 



1 k J^_A 



h/[om/£; + a] 



n=0 m: 



^ m/k + ak r 



=k~ n 



-EE 

n=0 j=l 



1 



j + afc n 



In' (bj + a) 1 g ln € [&(./' - m/A-) + «] 



m=l 



j + m/k + ak T 



(1.96) 



&=fc- 



Our results have application to Dirichlet L functions, as these may be written 
as a linear combination of Hurwitz zeta functions. For instance, for \ a principal 
(nonprincipal) character modulo m and Re s > 1 (Re s > 0) we have 



X(k) 



1 



k 



(1.10) 



fc=l 



fc=l 
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Proof of Propositions 

Proposition 1. We apply 
Lemma 1. We have for Re s > 1, 



S 1 oo oo 



C(s,a) 



a- s 1 



s- 1 2 4 



n=0 j=0 



+ 



1 



-r oo oo 

= a- s + - £ 2^- 1 ) £ 
2 n= o j=o 

We start from the representation 



(2~ n j + a) s (2~ n j + 2~ n + a) s (2~ n (j + 1/2) + a) £ 

1 



(j + 1 + a2 n ) s (j + 1/2 + a2 n Y 



. i s a a 
C(s, a) = — + 



1-s 



-<ix, Re s > 0, 



(2.2) 



2 s - 1 Jo (x + a) s+1 
and employ the functions [5] /(x) = — Pi (a;) and g%(x) = f(x) — f(2x)/2, with 
En=o92(2 n x)/2 n = f(x). We find 

g 2 (2 n x) 



C(s,a) 



a a 



a r 



(x + a) ; 



+i 



<ix 



a a 
,i 



l n=0 -— n 



,tS^ (2-"y + a) 



s+1 



dy 



j+i/2 (2- n y + a) 



s+l 



(2-3) 



Carrying out the integrations then gives the first line of (2.1). Elementary manipu- 
lations then yield 



C(s,a) 



a~ s a 1 ' 3 1 00 



2 



n=0 



1 1 

+ 



a~ s a 1 - 8 1 00 



+ + - V 2 n(s - 1} 

s - 1 4 f< 



1 



(j + l + a2 n ) s (j + a2 n ) s (j + l/2 + a2 n Y 
1 1 



n=0 



a s 2 r 



2E 

j=0 



(j + l + a2 n ) s (j + l/2 + a2 r - 



(2.4) 
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giving the rest of the Lemma. We remark that (2.3) is generalized in (2.14) for values 
of k > 2. 

We now use standard expansions about s — 1, including 



( j + P + o2 



n\s-l+l 



= (j + P + «2 n ) exp[(s - 1) ln(j + /3 + a2 n )\, (2.5) 



to write 

((s , „) - » = f) t*£ ln< a(s - !)«-! + 1 f; W a(s - 1)' 

1 ~ « ~ (s - If j [n In 2 - ln(j + 1 + a2 n )f [n In 2 - ln(j + 1/2 + a2 n ) 
^ o 2^ 



2 ^ ^ ^ P\ 

^ n=0 j=0 £=0 



j + 1 + a2 n 



j + l/2 + a2 n 



Comparing to the expansion (1.1) we have 

le (a) = - TTT + -In a 
+ (-1)* g g | [ n In 2 - ln(j + 1 + a2 n )]^ [n In 2 - ln(j + 1/2 + a2 n )]< 



n=o,=ok J + l + a2' 

ln m a 1 



7Tr + -Wa + EE 

^ 1 U n=0 j=0 



lnMa + ^1 



j + 1/2 + a2 
In* (a 



2 n 



ln £+1 n 1 oo oo 

T 1 " 71=1 j = l 

This completes the Proposition. 

We can rearrange (2.1) in the form 

From this follows Corollary 1. 



2j + 2 + a2 n+1 2j + 1 + a2 n + 1 



ln^a+^r) ln £ (a + ^ 



2j + a2 n 



2j-l + a2 n 



(2j + a n ) s (2j - 1 + a n Y 



(2.6) 



(2.7) 



(2.8) 
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Remark. It is easily verified that from Corollary 1 we may obtain the property 
d a ((s,a) = —s((s + l,a). Likewise it follows that ^(s, 1/2) = (2 s — l)C(s), and we 
provide a demonstration of this property from (1.5). We have 



1 \ os — 1 oo oo / i \j 

a s ,-)=2 8 + - — + y2"^v , 1 ' - 



(2.9) 



We transform the latter double sum to 



2(n+i)(s-i) ( z¥- = 2 S_1 

^ ~ j ; (.1 ■ 2'T 



n=0 



(-1) 3 ' f 2 n (s -i) f 
,t!0' + l) s £i ,tt(j + 2")^ 



= 2 



s-l 



c( s )-i-2 l - s c( s ) + E2 n(s ~ 1) E 



71=1 



i u + z n y 



(2.10) 



Then from (2.9) and (1.5) with a = 1, 



c, ( *, - J = 2 s + — T + r-\(s) - 2- 1 - C(a) + 2- 1 E 2" (s - 1} E 



IS — 1 OS— 1 



= 2 s + 



s-l 
= (2 s -l)C(s). 



+ 2 S C( S )-2 S -C(s) 



(2.11) 



Above, we used the alternating form of the zeta function so that 



-IE- 



-iy 



1 = (1 - 2 1 -)C(a) - 1. 



(2.12) 



Corollary 2. The first expression for 7 follows by performing one of the sums in 
the Brun result (1.6), 



7 = i-/3(i) = i-EEVt7 

n=lj=l Z 1 "J 



(2.13) 



In order to obtain the second expression for 7, we apply the duplication formula of 
the digamma function, 2ip(2x) = 2 In 2 + ip(x) + ip(x + 1/2). 

Remark. The second expression for 7 in Corollary 2 is precisely (2.24) of [5] with 
a = 1. 

Corollary 3. We have from (1.3) 

1 1 



7o (a) = -V>(a) = 
Summing over j then gives 



1 1 00 

In a H h - > 

a 2 ^ 



j + l + a2 n j + l/2 + a2 n 



(2.14) 



1 1 00 r / i\ 

70(a) = -lna + - + - £ ^a2 n + - -^ (a2 n + 1) 

a ^ n=0 L ^ ^' 



(2.15) 



We then use both the functional equation ip(x + 1) = ip(x) + l/x and the duplication 
formula of the digamma function to obtain 

1 



^ 00 

70(a) = - In a H 

a ^0 



i> (a2 n+l ) -^(a2 n ) - In 2 



a2™ +1 



(2.16) 



The Corollary follows. 

Proposition 2. We simply outline the proof. We start again from the integral 
representation (2.2). Now we make use of the functions for k > 2 and f(x) = —P\{x) 



0, 



9k{x) = f(x) - -f{kx), 



(2.17) 



with En= ^ = /(*)• Then 



ct s a} s 
C{S > a) - — - — 



V — 



n=0 



g k (k n x) 



k n Jo (x + a) s+1 



-dx 



y — 



n- 



9k{y)dy 



" Q k 2n Jo (k- n y + a) s+1 



™J_™ri+l g k (y)dy 
S 2^ L2n 2^ /, 



OO oo 

= S E T^Tt E 
n=0 ^ j=0 



If IN /-J+V* If 3\ f 
2V-k)S, + 2V-k)l, 

dy 



2 \ fc/ /j+i/fc 2 V fc/ Jj+2/k 



1/1 \ P+ 1 

2 U / Jj+(k-i)/k 



{k~ n y + a) s+r 



(2.18) 



In the last step we have used the values of g^{x) on sub intervals 
1, 2, . . . , k, 

. . 1 / IN (7-1) 



i - 1 i 



(2.19) 



fc ' k) ' 

In particular, the difference of these values on consecutive subintervals is simply 1/k. 
We then perform the integrations and collect the terms to find for Re s > 1 



a s a 1 s 



00 1 °° ( 1 / 1 
n =o ^™ j=o I 2 VA; 



C(s ' a) = — + s -l 



+ 



fc-i 



(2.20) 



(6j + a) s (6(j + 1) + a) 

Expanding about s = 1 and using the definition (1.1) gives the Proposition. 

Summary 

Among other results, we have found computationally useful series representations 
of the Stieltjes constants. Although the series representations are doubly infinite, 
one of the series converges exponentially quickly with parameter k > 2. We have 
obtained in (2.16) a parameterized representation of the Hurwitz zeta function. In 



9 



the very special case of k = 2 and a = 1, we recover the much earlier result of Brun 
[1] of a double series expression for the Riemann zeta function. Since 70(a) = —ip(a), 
where ip is the digamma function, even at the lowest order, we effectively have found 
series representations of the harmonic numbers, generalized harmonic numbers, and 
other mathematical constants. We recall that the generalized harmonic numbers 
= Y?k=i l/^ r may be readily found from the polygamma functions ip^' 1 ^. Very 
special limit cases of our results include parameterized series representations for the 
expressions £(0, a) = 1/2 — a and ('(0, a) = lnT(a) — (1/2) ln(27r). Further, our results 
have application to series representation and expansions of Dirichlet L functions. 
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